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Abstract — The necessity of accurate channel estimation for 
Successive and Parallel Interference Cancellation is well known. 
Iterative channel estimation and channel decoding (for instance 
by means of the Expectation-Maximization algorithm) is partic- 
ularly important for these multiuser detection schemes in the 
presence of time varying channels, where a high density of 
pilots is necessary to track the channel. This paper designs a 
method to analytically derive a weighting factor a, necessary 
to improve the efficiency of interference cancellation in the 
presence of poor channel estimates. Moreover, this weighting 
factor effectively mitigates the presence of incorrect decisions at 
the output of the channel decoder. The analysis provides insight 
into the properties of such interference cancellation scheme and 
the proposed approach significantly increases the effectiveness 
of Successive Interference Cancellation under the presence of 
channel estimation errors, which leads to gains of up to 3 dB. 

I. Introduction 

Multi User Detection dMUDb has made very signifi- 
cant steps from conceptual tool into practice [l]-[4]. Itera- 
tive Channel Estimation ( ICEb and Interference Cancellation 
(IC) [5], [6] (for instance by means of the Expectation- 
Maximization algorithm [7]) is one of the most appealing 
techniques, partly due to its linear complexity in the number of 
users. According to this principle, every time new estimates 
about the nansmitted symbols are available, the channel is 
estimated again so as to refine the precision of the Channel 
State Information (CSI). On the other hand, as new (hopefully 
more accurate) CSI is available, the MUD is repeated, lever- 
aging the CSI improved quality. Such process is iterated until 
convergence is achieved. In addition, the performance of IC is 
enhanced by the usage of soft-estimates of the coded symbols 
based on the computed symbol A Posteriori Probabilities 
( lAPPk ). If the bits are not correctly decoded, the estimated 
bits are regarded as unreliable and hence little or no signal 
will be subtracted. 

This idea has proved to be very effective for slowly varying 
channels, because reliable channel estimation is possible [5], 
[7], [8]. The usage of data symbols for the channel estimation 
process is in principle even more atnactive for time varying 
channels, because this approach increases the sampling fre- 
quency of the channel and enables better tracking thereof. 
Such type of channels shows up for example with mobile 
users or when phase noise is non-negligible. The impact of 
phase noise is particularly relevant for instance for applications 
at high carrier frequency (say Ku or Ka band, typical in 
modern satellite environments), where high stability oscillators 
can be very expensive; thus cheap, consumer-grade terminals 
are significantly affected by such problem. A pilot based 
CE would require in these conditions a too high density of 
known symbols to accurately track the channel. The problem 



of reliable channel estimation is even more pressing with MUD 
based on Successive- or Parallel- Interference Cancellation, 
since channel estimation errors lead to residual, non-cancelled 
interference. While a large number of pilots would entail an 
excessive overhead, the amount of symbols used for channel 
estimation purposes may be increased by using the output 
of the decoding process and therefore iterative CE-MUD 
approaches sound very suitable for this task. 

On the other hand, a subtle issue with iterative channel 
estimation MUD seeps in. The goal of channel decoding is 
to find the sequence of coded bits that best fits the received 
samples, given a set of known and fixed channel estimates. If 
the process of channel estimation and data detection becomes 
iterative, the receiver will find the sequence of symbols 
and channel estimates that jointly fit the received samples. 
Therefore, the number of degrees of freedom to "explain 
the observations" increases and hence the decoded bits will 
have large lAPPfe and will be regarded as reliable even when 
they do not correspond to the actually transmitted bits. The 
effectiveness of soft interference cancellation will be severely 
limited by this problem. 

Section [TT] will present the system model, while the main 
contribution of this paper is further described in Section [Til] 
and consists in the definition of a modified soft interference 
cancellation scheme that is especially effective for time- 
varying channels. The key idea is to multiply by a coefficient 
a the reconstructed waveform to be subtracted. Successive 
Interference Cancellation (ISICt is in fact an iterative algorithm 
to find the solution of a linear equation system. Therefore, it is 
equivalent to the Gauss-Seidel iteration. Applying a weighting 
factor is equivalent to Gauss-Seidel over relaxation, and it is 
known to improve convergence [9]. The idea of applying a 
weighting factor is not new per se, as it has appeared in [10], 
[11]. However, the derivation of the value for this IC factor a 
was heuristic and driven by numerical simulation [11]. Instead, 
an analytic derivation of such factor is provided and hence it 
is possible to optimize its effectiveness and get insight into 
the properties of this IC technique. Moreover, our proposed 
approach works also for asynchronous users and time varying 
channels (which is not the case of [7]). Section[IV]will evaluate 
the performance of this scheme by means of simulations and 
finally Section W\ will draw the conclusions. 

II. System Model 

The first part of this section will introduce the essential 
paper notation. Then, the basic model that will be further 
analyzed in Sections [TTtJ and [TV] is described. The final part 
(Section III-Cb will outline some important extensions to the 
basic system. 
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Fig. 1. The architecture of a generic stage of the multistage SIC receiver. 



A. Essential Notation 

Let us consider a generic complex number a. Its real part, 
magnitude, argument and conjugate are denoted as Re(a), \a\, 
Arg(a) and a*, respectively. The estimate of any quantity a 
(complex or real) is denoted as a. The Hadamard product (i.e., 
elementwise) between two vectors a\ and is denoted by 
a\ ■ a.2- Finally, the expectation of a random variable a is 
denoted by E[a\. 

B. Basic Model 

The system is composed by U users that transmit simul- 
taneously and in the same frequency to a common receiver, 
all involved nodes are equipped with a single antenna and 
no spreading is used (a possible extension to multiple an- 
tennas will be outlined in the next subsection). Each user 
u £ {0, . . . , U — 1} generates a sequence b u of B u inde- 
pendent, identically distributed information bits. Each set of 
information bits is channel encoded into C u coded bits by 
means of a rate R u = B u /C u channel code. In the simu- 
lation section the UMTS Turbo code will be assumed [12], 
but in fact the proposed algorithm does not depend on the 
specific code, as long as a Soft-Input Soft-Output dSlSOb 
decoder is adopted [13]. The coded bits are interleaved and 
then modulated by means of an M-PSK modulation into 
Xd,u = C«/ ]og 2 (M) modulated data symbols, which are time 
multiplexed with X p<u pilots symbols for channel estimation 
purposes. Let us denote by X u the whole set of data and pilot 
symbols. They are pulse modulated and transmitted through 
a frequency flat AWGN channel, which adds complex zero 
mean circularly symmetric White Gaussian Noise (WGN) with 
per dimension variance a\. The channel is not assumed to be 
slowly time varying, i.e., it may change from symbol to symbol 
due to, for example, phase noise or Doppler. We shall assume 
that the former is indeed present due to the instability of the 
transmitters and receiver oscillators. 

The signals from the U different transmitters arrive at 
the receiver. The users will be assumed in the description 
of the system model to be symbol and frame synchronous 
for notational simplicity. This assumption will also make the 
derivation and discussion of the a factor easier to follow in 
Section |III1 However, all simulations have been carried out for 



the symbol and frame asynchronous case, thus the performance 
evaluation is realistic in this respect. 

The receiver performs matched filtering and sampling and 
obtains a sequence of samples y t , where t is the discrete 
temporal index, < t < m&x({X u }) — 1: 
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where h t _ u = |ft.f, iM |e l6lM ' is the frequency flat complex channel 
at time t for user u, Xt >u is tne modulated symbols transmitted 
at time t by user u and w t is the complex WGN. In our 
case the variability of the channel coefficient ht, u is induced 
by the phase noise due to oscillator instabilities. However, 
the analysis presented is valid for any flat complex channel. 
A random-walk (Wiener) phase noise model with Gaussian 
increments is adopted [14]. The additional phase due to such 
impairment can be described as: 
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The increment A t u is a zero mean real Gaussian random 
variable with variance at. The receiver performs a multistage 
ISICI with iterative CE and Channel Decoding (CD), which is 
depicted in Fig. Q] In order to improve the error correction 
performance, the users are ordered in descending order of 
power [1] and we assume without loss of generality that the 
user index u corresponds to the decoding order. At each stage 
s, < s < 5 — 1 the u-th user is decoded based on the 
following signal: 
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which is equal to the received samples minus the estimated 
interference from the other users. The interference from the 
other users is reconstructed using the A Posteriori Probability 
(lAPPt rather than extrinsic probabilities. It has been shown 
that when a graph based approach is used for IMUDl extrinsic 
probabilities should be used rather than lAPPl [151. However it 
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is not known whether [APPI or extrinsic probabilities should be 
used in case a different IMUDI algorithm is used, such as ISICI 
[16]. In our case lAPPl provide better performance. Note that 
Vt,o,o = Vt- The estimates of the channel and of the modulated 
symbols from stage s of user u at time t are denoted by h t . SyU 
and Xt lS ,u- Estimates from the same stage s are employed 
for users 0, ...,u — 1 already decoded in this stage, whereas 
estimates from the previous stage s — 1 are used for the other 
users u + 1, U — 1. Note the presence of the factors a s . u , 
which will be the focus of the next section. Their purpose is 
to weight the confidence on the estimated waveform. 

The values for ht lS ,u and Xt s ,u are derived from the 
channel estimation and channel decoding, respectively. In 
order to fight the channel variability, an iterative channel 
estimation/channel decoding algorithm which follows the 
Expectation-Maximization (EM) approach is adopted [7], 
[17]— [19]. The derivation of the steps is carried out according 
to [18] and it can be shown that the E-step is simply the 
channel decoding. From the E-step, the lAPPfe of the coded 
bits ct lU can be evaluated [13] and hence E[a?t )Sl „] Q and 
E[\x t , s .u\ 2 ] are computed for the sake of channel estimation 
and interference cancellation. Moreover, the channel decoder 
gives as additional output Xt, a ,u> which can be (but need not 
be) E [xt iS . u }, i.e., the Minimum Mean Square Error ( IMMSEI i 
estimate of Xt. SjU . In addition, E [|xt, s ,u| 2 ] will be denoted 
with a slight abuse of notation |£t,s,«| 2 - Note that also known 
symbols (for instance the pilots) are employed in the CE 
process. Since these symbols are known with no uncertainty, 
E[xi, s ,„] = x t , u and E[|x 2 iS> J] = |xf jU |. 

On the other hand, the M-step corresponds to the channel 
estimation carried out according to these equations [18]: 
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Note that the channel is time variant and therefore all 
estimates are performed over a sliding window of size 2W+1 
samples around the desired time index t, where W depends 
on the coherence time of the channel. 

After the S MUD stages, the channel decoding estimates 
of the information bits b u are hard thresholded and given as 
output of the MUD/CD process. 

C. Extensions 

We briefly outline here two extensions of this framework. 
The first one concerns symbol asynchronous users and the 
latter V-BLAST [20]. 

In the first case, the users are symbol asynchronous and 
the receiver first samples the overall signal into discrete time 
samples {y t } with a timing not necessarily related to the 

'The expectation is actually conditioned to the observation, E[xi >s>u |j/t]. 
In the following the conditioning on the observation will be dropped for 
notational simplicity. 



sampling time of the users. Such timing will be called "the 
common time frame" and is opposed to the optimal matched 
filter sampling choice for each user, that will be called the 
"optimal single user time frame". The channel estimation and 
decoding block receives samples at the common time frame, 
which are interpolated to the optimal single user frame. All 
operations that are related to the tagged user (e.g., CE and 
CD) are performed in the latter system and the computed 
channel and symbol estimates are then interpolated back to 
the common time frame. After this operation, the reconstructed 
waveform can be subtracted from the received samples y t . 

The second extension concerns V-BLAST. The receiver and 
the transmitters are assumed to be equipped with N r and 
M t antennas, respectively, and the transmitters send inde- 
pendent streams. The users are ordered according to some 
usual criterion (say, maximum post-processing SNR) and the 
received samples are first processed by means of Zero Forcing 
or IMMSEI filtering and then the selected user is decoded by 
means of the previously described iterative channel estimation 
and decoding approach. Once the estimates are ready, the 
weight factor a s>u is computed, the waveform of the detected 
user is reconstructed, multiplied by a S;U and subtracted from 
the received samples. This approach is reminiscent of Turbo 
MIMO [21], with the difference that also channel estimation 
is iterative and for the presence of the a factors. 

III. Optimal Partial Interference Cancellation 

Let us denote as y u = h u ■ x u the noiseless signal received 
from user u, being x u and h u the symbols and channel 
coefficients from user u. Stage and temporal indexes are 
dropped here for notational simplicity. Let y u = h u ■ x u be 
an estimate of y u . Assume now that we want to cancel the 
interference caused by user u on other users. In accordance 
to the model defined in Section [EI] a ISICI scheme is applied in 
which the estimates are first multiplied by a weighting factor 
a u and then subtracted, with < a u < 1. The interference 
cancellation efficiency [22] for such a scheme can be defined 
as: 

\Vu - ot u y u 
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In an ideal case, when the interference caused by y u is 
completely removed, (3 = 1. On the other hand, if no 
interference is removed, /3 = 0. 

If the adopted modulation is PSK, the interference cancel- 
lation efficiency can be expanded as: 
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The assumption of a PSK modulation implies that \x u \ 2 — 
l,Vx„. Extension of such analysis for non-constant envelope 
modulations is left for future work. 

The last step from Eq. (0 to Eq. (TTOb is justified because 
the last term of Eq. (0 is generally zero under the reasonable 
assumption that the channel estimation error is uncorrelated 
to the value of the channel. Under these conditions: 
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M 2 



is the normalized mean square 
channel estimation error. Mote that the second term in Eq. [10] 
is related to the reliability of the channel decoding, while 
the third term conveys information on the channel estimation 
accuracy. 

Setting the first derivative of f3 with respect to a u to 0, it 
can be easily verified that the value a Uy0p t that maximizes [3 

E[Re{x u -x* u }} 
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Note that the optimal value of the weighting coefficient 
ct u ,opt depends on E [|x u | 2 ] , E and E [Re{x u ■ x* }]. The first 
term can be computed at the receiver, since x u is available at 
the output of the ISISOI decoder. It can be assumed that the 
second term can be estimated given the SINR at which the re- 
ceiver operates. The computation of E [Re{x u ■ x* }] involves 
perfect knowledge of the transmitted symbols. However, a 
good estimate of this term can be obtained replacing x u by 
hard decisions on x u . E can be also reliably estimated by 
means of a look-up table based on an estimate of the SINR. 

Furthermore, note that the right hand side of Eq. (fT~2t can be 
described by two terms, the first depending on the reliability 
of the channel decoder output, while the second depends on 
the accuracy of the channel estimation. When pilot aided |CE| 
is considered, the optimal weighting coefficients for data and 
pilot symbols are different. In case of data symbols a op t is 
given by: 
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where xffl are the data symbols of user u. In case of the pilot 
symbols given that their value is known a priori, Eq. (fT3l can 
be simplified to: 

a { P ] =-^ (14) 
1 + E 

Simulations verified that al„ pt > ai„ pt , with equality for 
high SNR. For low SNR, the value of the data symbols can 
not be reliably estimated , x u and Xu tire uncorrelated and 

0. On the other hand, when 

(d) , 



therefore: 
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the estimates become more reliable, E 

E [|xi d) | 2 ], and therefore a { ^ pt ~ a$ pt . 

^et us now consider the case in which x u = x u . This is 
the case of pilot symbols, or data symbols when we operate 
in a region where the Bit Error Rate (BER) ~ 0. Let us now 
assume that the normalized channel estimation error, E, is 
known. We shall call full ISICI the conventional ISICI scheme 



with a = 1. Under full [SIC] the normalized interference caused 
by user u to other users can be expressed as: 

"f itr 1 2/u — 2/u I w | h u ■ x u — h u ■ x u - 
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(15) 

In a system with optimal partial ISICI according to Eq. (fl2l) . 
the normalized interference can be expressed as: 
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We can define 7 as the ratio of the normalized interference 
with full ISICI and optimal partial ISICI 
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Note that since E > 0, 7 is a monotonically decreasing 
function of E and 7 < 1. 

Let us now consider the opposite case, where the channel 
estimates are perfect, but the estimates of the data symbols 
are not. Let X u be the Log-Likelihood Ratios ( ILLRfe ) of the 
lAPPfe obtained at the output of the ISISOI channel decoder of 
user u. For BPSK, the symbol estimates are calculated as: 



tanh 



(18) 



This choice is optimal from a IMMSEl consideration [5]. In 
other words, under perfect Channel State Information ( ICSIb 
and using soft decision according to Eq. ( fT8l , full ISICI (a = 
1) minimizes (3. This was verified by means of simulations. 
However, when Expectation-Maximization (IEMI ) based ICEI is 
used, this is no longer true. Fig. [2] shows the evolution of the 



BER vs. E 



.(d) |2 



for perfect ICSJ and EM based ICEl using 
the UMTS turbo code with code rate 0.53. The figure shows 
results for a single user case and a two user case according 
to the setup described in Section [TV] (user 2). For the single 
user case it can be seen how for the same BER, E lii^! 2 
is higher for EM based ICEI than for perfect ICSIl This means 
that the system is overconfident in the decisions it takes. The 
figure also shows how in the presence of two users E \x u d ^ | 2 
is higher than in the single user case, both with perfect ICSIl 
Finally the two user case with EM based ICEI produces the 
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In Section [TV] it will be shown how optimal partial [SICl can 
translate into a considerable performance improvement in a 
system where ISICI with IEMIICEI is used. Especially when the 
power imbalance between users is low and hence the channel 
estimation error is high, full ISICI is unable to bootstrap due to 
the high residual interference from other users. Optimal partial 
ISICI decreases the residual interference enabling gains of up 
to 3 dB for power imbalance 1 dB. 

IV. Simulation Study 

A simulation study was carried out in order to evaluate the 
performance of the proposed partial interference cancellation 
scheme. A simulation setting with two symbol and frame 
asynchronous users is used. The channel code is the UMTS 
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Fig. 2. BER vs. soft values in a 1 user system (solid lines) and 2 user system 
(dashed lines) with 2 dB power imbalance. 

turbo code [12]. The code rates of user 1 and user 2 are 0.72 
and 0.53 respectively^ The number of ISICI stages is set to 7 
in all cases, as no meaningful performance improvement is 
observed after this point. In each stage each user performs 
15 iterations of the IEMI algorithm. Both transmitters employ 
the same random interleaver. The modulation used is BPSK 
for both users, and the code block length is C u = 5000 
symbols, Vit. The number of pilot symbols per code block 
is X UtP = 256, Vu. Pilot symbols are boosted by 3 dB with 
respect to data symbols. The channel is assumed to be AWGN 
and the phase noise increase between two consecutive symbols 
has a variance = 100//o rad 2 /symbol, where fa is the baud 
rate (in our simulations 10 kBaud)0 Different simulations were 
carried out for power imbalance, Pjjui u 2 > ranging from 1 dB 
to 7 dB, being user 1 always the strongest one. In contrast to 
much of the work on IMUDI no spreading is employed^ 

Our scheme is compared to multistage full ISICI with EM 
based CE (very akin to [7]). Fig. [3] and Fig. [4] show the BER 
of user 1 and 2, respectively, for a power imbalance of 2 dB. 
If we analyze the SNR at which the different schemes reach 
a BER of le — 4, for both users optimal partial [SlCl results in 
a SNR gain of 2 dB with respect to full ISICI Moreover, we 
can see how optimal partial [SICl with 2 stages outperforms full 
ISICI with 7 stages. However, both ISICl schemes with estimated 
ICSII are far from the performance achieved by perfect ICSll For 
optimal partial ICSll the loss with respect to perfect CSI is 2 dB 
for both users, half of the gap of the state-of-the art system. 

Fig. [5] shows the evolution of E with the SNR for user 2. It 
can be seen how after 7 ISICI stages the scheme with optimal 
partial [SICl outperforms full ISICI However the performance of 
both schemes is far from the performance of the single user 
case. The reason for this is the remaining interference from 
user 1. It is also quite remarkable that optimal partial [SICl after 
2 stages outperforms again full ISICI with 7 stages. Another 

2 Such rates corresponds to the capacity of a BPSK constrained AWGN 
channel with nominal SNR for the second user 1.5 dB and power imbalance 
2 dB. The SIR of the second user is computed as 2 dB, that is to say, IC 
of the first user leads to a 5 dB improvement. Finally, the channel codes are 
assumed to be 1 dB away from capacity. 

3 This statistics has been derived from satellite transceivers normally avail- 
able to the authors. 

4 In spite of the use of a low cardinality constellation like BPSK and the 
fact that no spreading is used, the interference caused among users can be 
considered Gaussian due to phase noise and the effect of the interleaver. 
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Fig. 3. BER vs. SNR for the first user, power imbalance 2 dB. 
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Fig. 4. BER vs. SNR for the second user, power imbalance 2 dB. 

important fact is that there is a floor at E = 0.031, which is 
caused by the presence of strong phase noise. Such problem 
clearly limits the performance of ISICI since the interference 
can never be completely canceled. The evolution of E for 
user 1, which is not shown in any figure, is very similar to 
that of user 2, presenting also the same floor. 

The evolution of ot^P with respect to the SNR is shown in 
Fig. [6] It can be seen how a\f ' increases with the number of 
ISICI stages. The figure also shows how the maximum value of 
a4 is 0.97 < 1. The maximum possible value (i.e., 1) is not 
reacheable because of the presence of a floor in the channel 
estimation. The proposed scheme is aware of this floor and 
avoids full removal of the estimated signal. 

Finally, Fig. [7] shows the SNR at which optimal partial [SlCl 
and full lSICI reach BER le — 4 for different power imbalances. 
It can be seen how the gain provided by optimal partial ISICI 
is high for low power imbalance and decreases as the power 
imbalance increases. Note that in many cases the optimal 
choice of a halves the gap to perfect CSI given this time 
varying channel. 

V. Conclusions 

The principle of soft ISICI hinges on the idea that a low 
confidence on the decoded bits would reduce the amount of 
subtracted interference. The EM algorithm is known to be an 
effective tool to improve the performance of MUD. However, 
it has the side effect of leading to overconfidence in the ILLRk 
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Fig. 5. Normalized mean square error in the channel estimate E for the 
second user with a 2 dB power imbalance between the two users. 
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Fig. 6. Evolution of the a factor for the second user as a function of the 
SNR. The power imbalance is 2 dB. 



which has a negative effect on the performance of ISICI The 
optimal weighting factor a for ISICI has been derived, which 
is able to partially compensate for this overconfidence in the 
ILLRfe . Moreover, it is also able to mitigate the negative impact 
of imperfect ICSII The proposed technique shows gains up to 
3 dB when the detected users are received with comparable 
powers. Furthermore, the performance gap between perfect 
CSI and the state of the art is roughly halved by the suggested 
scheme. 
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